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ABSTRACT 

Formulae which have been proposed for representing the longi- 
tudinal profiles of cosmic ray air showers are compared, and 
the physical interpretation of the parameters they contain is 
examined. Applications to the problem of energy calibration 
are pointed out. Adoption of a certain especially simple for- 
mula is recommended, and its use is illustrated by means of 
examples. 


1. Introduction. At all primary energies above ^ 10 3 GeV, incoming cos- 
mic ray nuclei deposit more than half of their energy via the soft compo- 
nent. At the highest observed energies the fraction E E ^/E approaches 
90-95%, where E EM is the energy deposited by the composite electron-pho- 
ton cascade resulting primarily from tt 0 decay y-rays. The single most 
important step, therefore, in determining the cosmic ray energy scale for 
E > 10 6 GeV, where experiments above the earth* s atmosphere currently 
give way to experiments using air showers, is measurement of E EM . 

The use of simple mathematical formulae to represent the lateral 
structure of extensive air showers began at an early stage in the study 
of shower phenomena (cf . Bethe quoted by Williams 1948) , and such formu- 
lae continue to be useful tools. Analogous formulae for the longitudinal 
structure developed much later, possibly because accurate experimental 
data have taken much longer to obtain, and possibly because analytical 
methods have tended to be supplanted in recent years by Monte Carlo cal- 
culations. 

With the advent of techniques for directly measuring the longitudin- 
al profiles of individual air showers (Hammond et al. 1978, Grigoriev et 
al. 1979, Cassiday 1981, Cady et al • 1983) it seems an appropriate time 
to re-examine formulae that have been used to describe these profiles 
(Greisen 1956, Linsley 1967, Longo and Sestili 1975, Gaisser 1976 and 
1979, Gaisser and Hillas 1977, Sass and Spiro 1978, Dyakonov et al. 1981). 
By analogy to formulae used for describing the lateral structure of vari- 
ous air shower components I will call them 'trial functions' . In case of 
lateral structure such functions have had an important role in comparing 
results of different experiments. One expects longitudinal trial func- 
tions to be useful in the same way. In case of the lateral structure of 
electrons, trial functions are also used in the conduct of experiments, 
to derive from raw data on shower density a global measure, the shower 
size N, as well as core location, age, and measures such as S^qq* The 
corresponding step in dealing with profile data is to find N m , long re- 
cognized as being one of the most reliable estimators of primary energy 
(Clark 1962), from data on N vs x. This step also yields x m , the shower 
elongation, and may also yield the profile width cr x . In case of the 
electronic lateral structure, . trial functions are formulated in such a 
manner as to reflect the physical processes that govern cascade develop- 
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ment. Thus core distances are expressed in Moliere units, and steepness 
is controlled by a parameter relating plausibly to shower age. One 
aspires to achieve similar success in case of the longitudinal structure. 

2. Electromagnetic cascades. The longitudinal trial function proposed by 
Greisen (1956) provides a good illustration: 

N = (0.31/6 Q ^)exp[t(l - — Ins)], where s = 3t/ (t + 2$ q ) . (1) 

Clearly its structure was influenced in several ways by theoretical re- 
sults obtained with diffusion equations (see Rossi 1952) : 

1) Thickness is measured in units of the radiation length (t = x/x Q ) . 

2) Primary energy is measured in units of the critical enerqy E (ft = 

ln(E/E c )). c ° 

3) There is a built-in elongation-energy relation, x = x 0 3 0 , which is 
approximately correct. 

4) There is a built-in N m relation, N m = (0.3lE/E c )/ [In (E/E c ) ] * , which is 
also approximately correct. 

5) The shower age s is incorporated in an approximately correct manner. 

The price^paid for (5) is mathematical inconvenience. Expressions for N, 
Nt and Nt cannot be integrated in closed form, so simple formulae for 
the track length integral, the average depth <x> and the profile width 
o x , provided by the exact theory, must be patched on ad hoc. 

An alternative form is the gamma distribution, recommended in the 
current Particle Properties Data Booklet for representing results from 
EGS, a well known Monte Carlo program for simulating electronic cascades. 
Gamma distributions have been chosen independently by several workers to 
represent air shower profiles as well. Indeed, the choice has been unan- 
imous, and this is the choice recommended here. 


3. Mathematical properties of gamma distributions. I will write this 
distribution in the form 

N = N £ g e~ q? , 
o ' 

calling the numerical constant q the 'index'. The optimum values of £ 
and N are given by * _ a 

E m ' 1 • l3 > H „ = V - 

while the mean value and variance are given by 

<5>= ii_t2 (5 , 

The normalization is given by /Nd£ = r(q+l)/q q+1 , 


( 2 ) 

(4) 

( 6 ) 

(7) 


which admits approximation using Stirling's theorem 

/Nd£ a N o (27T/q) i e“ q = N m (2i:/q) i . (8) 

I propose to determine the parameters of (2) by using 

1) the elongation-energy relation, x m vs InE (note from (3) that £=x/x ), 

2) a relation due to Kraushaar (1957) between the typical profile of an 
individual proton shower and the average profile of proton showers, 

3) experimental data on a x , and on E EM /E, the fraction of primary energy 
given to electrons. 
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First, however, I will show that (2) provides an excellent fit to elec- 
tromagnetic cascade profiles. 


4. Electromagnetic cascades revisited. Results of cascade theory in ap- 
proximation B are (Rossi 1952) : 

X m = X Q (1.01(3 + A), <X> = X Q (1.013 + B) , rr 2 = M Cl O ^ 


X Q (1.613 + C) , (9) 


where 3 = ln(E/E c ) and A,B,C are constants a, l whose exact values depend 
on whether the primary particle is an electron or photon. Using (3) and 
substituting in (6) I solve for the index 

q = 0.63 + D , (10) 

^fr®^® D — —0.4, 0.2 for primary electrons and photons, respectively. 
Substituting this result in (5) , I find 

<X> " X m = 1,6x o ' (11) 

expressing the fact that the profile is unsymmetrical . To conserve ener- 
gy the track length integral must equal E. Thus 

/(E c /x Q )Ndx = E . (12) 

Using (8) , which I solve for N , and again substituting for q and x 

m i m 

N m = (0.31E/E c )/[ln(E/E c ) - F] * , ( 13 ) 

where F = 1.7, 1.0 for primary electrons and photons, respectively. Re- 
sults (11) and (13) are almost exactly the same as given by the detailed 
theory. 

5. Applicatio n to air showers. One begins as in the preceding example, 
by letting 5 = x/x^, where Xjj^ is given by the elongation-energy relation. 
A reasonable choice for this is 

x m = A + D io lo 9 E (14) 

with A = 159 g/cm 2 and D 1() (elongation rate per decade) = 65 g/cm 2 
(Linsley and Watson 1981; E is in GeV) . The index q can be found as in 
the preceding example, using (6) , but first I will show an alternative 
method illustrating a deep connection between gamma distributions and 
cascade processes. 

It was pointed out by Kraushaar (1957) that the average number of 
electrons N at atmospheric depth x arising from a primary particle inci- 
dent upon the top of the atmosphere is related as follows to the average 
number of electrons at thickness x 1 below a nuclear interaction of 
such a particle: 

N(x) = f* N 1 (x')exp[-(x-x')/A)dx'/l , (15) 


N x (x) = [1 + A (9/9x) ] N (x) , ( 1 6 ) 

where A is the interaction mean free path of the particle, assuming that 
there are no fluctuations other than those in the starting level. Apply- 
ing these, it is readily shown that the condition for N and N]_ to be self 
similar} meaning in this case for one to be a gamma distribution if the 
other one is, is the followings 

* = x m /X * f 1 ?) 

Substituting (17) in (6) one obtains a very interesting result, 

a x = * 2 + Xx m / (18) 

which can be combined with (14 ) , giving 

- d 2 + AA) + AD 10 logE . 


(19) 
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This says that if the energy dependence of A is neglected, the energy de- 
pendence of has the same form as that of x^ just as it does in case 
of electromagnetic cascades. The energy dependence is described by a 
rate, analogous to the elongation rate, which is in fact proportional to 
the elongation rate. This new rate, which might be called the 1 width 
rate* or ’spreading rate ' , is equal to AD 1() . 

An experimental result on a x has recently been reported by the Fly^ 

Eye group. The value given, 220±33 g/cm 2 at 10 9 GeV (Baltrusaitis et al . 

1985) , agrees very well with a value derived from an average profile pub- 
lished by Grigoriev et al • 1983 for a slightly lower energy. Substitut- 
ing the average of these results, and values of A and D 1Q , in (19), one 
finds A = 58 g/cm 2 . 

This is appreciably greater than the mean free path for proton-air 
inelastic collisions found by other methods (see conference paper HEl.l 
-1) . It is reasonable to assume that A is greater than ^p a ,inel. for the 
same reason that holds in similar cases: neglect of development fluctua- 
tions. This needs further study. Pending results of such a study, a 
reasonable estimate of a 2 for energies other than 10 9 GeV is obtained by 
assuming that A is constant with the value found above. Then 

q = 3.36 + 1 . 121ogE , (20) and = 1.5-lo‘*+ 3.8*10 3 logE (21) 

The remaining parameter in (2) , N 0 / is found in the same manner as 
for electromagnetic cascades, using (7) or (8) , keeping in mind, of 
course, that the track length integral equals E^, not E. The relation 
between E EM and E is discussed in conference paper OG5.1-5 (or see 
Linsley 1983). A convenient formula representing results of that work is 

E em /E = 1 - 2 .8E“° * 1 7 (E in GeV) . (22) 
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